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Abstract 

The issue of quantum size effects of interactive electron-hole systems in spherical semiconductor 
quantum dots is put to question. A sharper theoretical approach is suggested based on a new 
pseudo-potential method. In this new setting, analytical computations can be performed in most 
intermediate steps lending stronger support to the adopted physical assumptions. The resulting 
numerical values for physical quantities are found to be much closer to the experimental values 
than those existing so far in the literature. 
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I. INTRODUCTION 



For about three decades, nanostructures — produced by many techniques such as etching, 
local inter-diffusion, particle suspension in dielectric media, or by self-assembly in matrices 
of a host material — are known as confined electronic systems because they are designed 
to restrict particle motion to a space region with a restricted number of dimensions. De- 
pending on their dimensionality, these structures are called quantum dots (OD), quantum 
wires (ID) or quantum wells (2D). The principal property of interest resides in their ad- 
justable quantized energy spectrum. The corresponding wave functions are localized within 
the confined region, but may extend over many periods of its boundary crystal lattice. A 
quantum dot (QD) is therefore a giant artificial atom which enjoys prospects for an increas- 
ing range of future applications: e.g. as a semiconductor laser as qubits for quantum 
information processing, as single-electron transistors in electronics, as artificial fluorophores 
for intra-operative detection of tumors, biological imaging or cell studies, etc. 

QD properties depend on their electronic structure. In fact, they contain a finite number 
of elementary charge carriers (from a few to a hundred), which may be conduction band 
electrons, valence band holes or excitons of the host substrate. There exists a vast number 
of aspects related to their own interactions as well as with the ambient electromagnetic field. 



Correlation effects 



^, phonon in- 



, interactions with boundary electromagnetic modes 
teractions j^, magnetic field effects 6], just to cite only a few topics, are currently vigorously 
investigated. 

In the early eighties, quantum size effects (QSE) occurred experimentally in spherical 
semiconductor QDs in the exploration of the optical properties of semiconductor microcrys- 
tals in an insulating matrix 7|, |8| . This phenomena has been observed in a large range of 
other confined structures: e.g. in quantum ribbons or in quantum disks 9|, in quantum 



wires [10| and, indeed, in quantum wells 



11| . It emerges in a widening of the semiconductor 



optical band gap. This is due to the increase of confinement energy for decreasing QD size. 
The valence band drops and the conduction band rises. Both effects constitute the leading 



contribution to the characteristic b 
strongly quantum- confined systems 



ue-shi 



12 



[ft in the semiconductors optical spectra of such 
13| . So, the confinement of oppositely charged car- 
riers does have significant effects on the electron-hole Coulomb interaction and, therefore, 
on exciton formation. 
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In order to apprehend correctly the origin of these quantum size effects, a first theoretical 
attempt to describe electronic properties dependent semiconductor QDs has been elaborated 
upon a particle-in- a- sphere model in the effective-mass approximation (EMA) 12, Is], [l^ . 
in which both electron and hole behave as non-interacting electrons and holes, trapped in 



a spherical infinite potential well but with different masses. Since 1973, the QD spherical 
shape has been used 1^ and continues to be very popular over the years 



0, Q, Q, Q, 



ITj. Because of the usual assumption of parabolic band structure, their effective mass is 



commonly defined through the inverse of the second derivative of their kinetic energy with 
respect to their momentum. Except for a pioneering work [3], the Coulomb interaction 
between electron and hole has been included and the excitonic contribution to the ground 
state energy has been taking into account by the Ritz variational principle. Some other 
authors have developed their own EMA model based on finite potential wells and improved 
agreement with experimental data for a significant range of QD sizes 



should also quote some other variational calculations, e.g. [20 



21 



Q, Q, ll9(. We 



I22I. In addition to 



spherical clusters, the case of cy 



and experimentally studied 18 



indrical shaped microcrystallites has been carefu 
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241 ] ■ as well as the case of quantum wires 25 



ly treated 



261 ]. Later 



on, more sophisticated models have been conceived. Some of them, called empirical tight- 
binding method (ETBM), have considered non-parabolic valence and/or conduction band(s) 
Others have used a reformulation of the so-called k-p perturbation theory including 



27 



non-parabolic bands [29 



30 



3l| based upon the Baldereschi-Lipari Hamiltonian [isl ]. 



However, despite the existence of numerous theoretical and/or empirical models, to the 
best of our knowledge, there exists actually no simple and comprehensive one, which offers 
a significant analytical treatment of the problem. This is the reason why this paper propose 
an alternative single EMA model for spherical semiconductor microcrystals, which allows 
us to clarify some tricky points. Moreover, we are able to exhibit the existence of an 
effective potential, whose presence significantly changes the QD ground state at small QD 
radius. Furthermore, in this formalism, it finally becomes possible, as an achievement, to 
analytically calculate some important quantities, such as e.g. the function rj introduced by 
Kayanuma in [l^, with a good agreement with experimental results. To this end, we shall 
introduce our model and recall some general properties in Sec. [HI Then, the two next Sees. 
mil and IIVI are devoted to the analysis of the two limiting regimes of semiconductor QDs 
which we shall explicitely describe and in which we carry out analytic computations in some 
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details. Finally, in Sec. |Vl the new effective potential mentioned above is introduced and 
we show how it leads to improve numerical results, much closer to experimental data. In 
the concluding section, we summarize our main results and indicate possible future research 
perspectives. 



II. QUANTUM DOT MODEL AND SOME GENERAL RESULTS 

Since our aim is to understand the semiconductor optical bands blue-shift as a main 
effect, we shall discard spin effects (electron- hole spin coupling or external applied magnetic 
field) and consider non-relativistic spinless electron or hole, trapped in a confining infinite 



spherical potential well. There exist other models with parabolic confinement ^, ll6|] or 
parabolic potential superimposed to an infinite potential well jl^l, which are used to explain 
certain spectroscopic data; but the concept of a QD size is then not so well defined. Here, 
we adopt the EMA model previously introduced. 



A. Intrinsic limitations 



First, the EMA assumption breaks down when a significant number of carriers wave 
functions inside the QD overlap. Thus, our results are expected to fail for very small 
nanocrystals. As matter of fact, at typical sizes of less than a hundred lattice spacings 
appear, in a semiconductor, magic numbers of clusters of which only some remain stable: 
e.g. nano crystalline sihcon stay solely coherent as clusters of Sii2, Siss, Si^g and 8145, if 



they contain less than 60 silicon atoms [3j]. When such size is reached, their band structure 
should be very deformed and, therefore, cannot satisfy the parabolic spectrum assumption 
on which the EMA is based. 

Second, by assuming spherical symmetry for the QD, we justify the splitting of the 
Schrodinger equation into a sum of a radial part and an angular part. This manipula- 
tion appreciably reduces the computation of the eigenstate wave functions and the energy 
eigenvalues. This approximation turns out often to be quite good since most synthesized 
nanocrystals possess an aspect ratio (defined as the ratio between the longest and shortest 
axes of the QD) smaller than 1.1. But, for high aspect ratio micro crystals, this spherical 
assumption becomes no longer valid. 
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Third, the real potential at the boundary of the QD is, of course, not infinite. Actually, 



there is a potential step with a standard height of magnitude from 1 to 3eV |2l|]. This value 
is generally quite large compared to typical electron and hole energies and, therefore, their 
tunnel conductivity from the nanocrystal to its surrounding should be neglected, except for 
very small QD sizes. Furthermore, the infinite potential well approximation implies that the 
charge carrier behavior inside a QD is totally insensitive to any externally applied potential 
or to the surrounding of the cluster. Although the surrounding effects may be sufficiently 
small to be neglected, the presence of a large external potential can in fact significantly 
modify the inside behavior of the microcrystallites. Thus, in order to test the validity of this 
approximation, we have to consider charged carriers isolated from the outside neighboring 
semiconductor. 



B. A free Quantum Dot model 



V{r) = V{r) 



Let V{r) be the confining potential well, defined in spherical coordinates as 

a 0<r <R, Region I; 
oo if r > R, Region II. 

Neglecting, for the moment, the electron- hole Coulomb interaction, in single parabolic band 
approximation, the Hamiltonian operator is (with ^ = 1) 

Ho = He + Hi, + Eg 



where m*j^ and ife,h denote the effective mass and the confined Hamiltonian respectively 
of the electron and of the hole, and Eg is the semiconductor energy band gap. The semi- 
conductor QD wave function can now be written as a product of its electronic and hole 
parts 

^'(re,rh) = 7/;(re)7/;(rh). 

The orthonormal eigenf unctions ipinm are labelled by three quantum numbers: / G N, n G 
N* = N - {0} and mG|-/,/]. 



R3l^ik^)y r "XR 
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where 

• YJ"(^, (p) is the spherical harmonic of orbital quantum number / and azimuthal quan- 
tum number m, 

• Jj,; (x) is the Bessel function of the first kind of index vi = l + ^ and variable x, 

• X[o,i?,[('^) is the characteristic function of radial interval, 

• {kin\in are the wave numbers in Region I, defined as the n**^ non-zero root of the Bessel 
function J;/;(x), thanks to the continuity condition at r = /?. 

The respective energy eigenvalues for electron and hole are expressed in terms of the same 
family of wave numbers {kin\i^ 

7-ie,h _ l^ln 

This clearly indicates that the continuum density of states of the semiconductor bulk should 
show atomic-like discrete energy levels with increasing energy separation as the QD radius 
decreases. 



C. Electron-hole pair Quantum Dot model 

From now on, we add the Coulomb interaction Vc between electron and hole to the 



Hamiltonian Hq, given by Eq. ([T]) and define the usua 



Hamiltonian H in the effective- mass approximation 13 . 



spherical shape semiconductor QD 
14,[l8|,Q as 



H = Hq + Vc{v,^) 

where k = Aire and e is the semiconductor dielectric constant. 

In order to simplify notations, we shall assume that Eg = 0. Deriving an exact analytical 
solution is arduous because of the Coulomb potential dependence in the electron-hole relative 
distance = |reh| = |rc — rh|, which explicitly breaks the spherical symmetry of the system. 
The common approach to this problem is to treat differently the interplay of the Coulomb 
interaction, which scales as (xR~^, and the quantum confinement, which scales as oc R~^. 
To handle these competing contributions, two regimes of electron-hole pair should be singled 
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out by comparing the order of magnitude of the QD radius R to the Bohr radius of the bulk 
Mott-Wannier exciton a* = -r— , fi being the reduced mass of the exciton. These are 

• the strong confinement regime, vahd for a size R < 2a* [l^, in which the potential 
well strongly affects the relative electron-hole motion, the exciton states consist then 
of uncorrelated electron and hole states; 

• the weak confinement regime, valid for a size R > 4a* [l^, in which the electron-hole 
relative motion and its binding energy are quasi left unchanged. The exciton could be 
treated as a confined quasi-particle of total mass M = m* + and its center-of-mass 
motion should be quantized. 

In any case, the Coulomb potential should be treated as a perturbation of the infinite 
confinement potential well. So, in order to evaluate the ground state energy of the exciton 
state using a variational procedure, we shall use the following wave function 

0(re,rh) = ^oio(re)V'oio(rh)0rci(reh), (3) 
where 0rci(rch) = <pTci{rch) = e"'^^"'^, a denotes a variational parameter, re,h = ke.hl, and 

V'oio(re,h) = V^oio(^e,h) = ^= sm — re,h • 

re,hV27ri? ^R ' 

The variational wave function ([3]) is the most natural choice we can make. The wave function 
V'oio(rc)^/'oio(rh) is simply the /ree Hamiltonian iJo ground state. Its presence in the trial 
function 0(rc, rh) insures the validity of the perturbation result to which the variational 
principle leads. Furthermore, in relative coordinates, the perturbation function f/'rcK'^ch) 
exhibits the exciton behavior, which is that of a hydrogen-like atom with a charge carrier 
mass \i. Because of this analogy, we expect the variational parameter a to depend on the 
Bohr radius as cr oc a*~^. Notably, in the weak confinement regime, the exciton ground state, 
which is a bound state, should mean that a ~ 2a*~^ ^ so that ^rcK'^ch) = e~^ — which is, 
up to a normalization constant, the ground state wave function of an hydrogen-like atom 
with its appropriate Bohr radius. 

Let us make a trivial remark which will be the basic argument for our later purpose: 
the spherical shape of the QD explicitly breaks the translation invariance of the Coulomb 
interaction. This remark on translation invariance and spherical symmetry breakdown, as a 
whole, suggests the use of a Fourier transform formalism in relative coordinates. Let T\f\ 



stands for the Fourier transform of the function /. Therefore, important matrix elements, 
such as the square of the variational wave function norm, the diagonal matrix element of the 
Coulomb electron-hole potential and the mean value of the electron-hole free Hamiltonian 
Hq in this variational state, may be advantageously computed with Fourier transforms of 
functions. 

The square of the norm of is then given by 

d^red\h^/'oio(rc)^/'oio(rh)0rci(rch) 



(27r)3' 

II — ^sin2(7rx)sin2(7ry) sinh(a/?x)e^'^^^ 
' (^JJv X y 

(4) 



where V = {{x,y) G M?/0 < x <y < 1}, while the different Fourier transforms are 



2 r^dx 



2 f dx 

•^[V'oiojlk) =J^J~ sin2(7rx) sm{kRx). 
Moreover, the diagonal matrix element of the electron-hole Coulomb potential is 



d^k 



kJ (2^)3-^Mk)^[^oio](k)^ 

= Z^1 II ^^sin2(7rx) sin2(7^^/)sinh(ai?x)e-"^^ 
kR aJJ-D X y 

(5) 

with v^rcifroh) = ^^^^^^ ^^^^ . Lastly, the mean value of the electron-hole Hamiltonian Hq can 

Teh 

be exactly evaluated as 



(010) 2/ii?2 gjj 2/ii?2 4 

^/ 

where the dimensionless variational parameter a' is defined by a = — and the exciton 

^ a* 

Rydberg energy by E* 



2/ia 



*2 ■ 



8 



As Eqs. (E]), (jlj), ([5]) and ([6]) are obtained without any further approximations — except 
the ones upon which the model was built — , they have to be valid everywhere whether it is, 
by construction, in the strong confinement regime or, by extension, in the weak confinement 
regime. Indeed, the Coulomb interaction itself should be always treated as a perturbation 
with respect to the infinite potential well, although its energetic contribution should not be 
maintained as a perturbation to the exciton confinement energy. Thus, the trial function 
global form ([3]) will remain acceptable in the weak confinement regime, but we will have to 
correct it by adding a further phase factor, which will depends only on the center-of-mass 
coordinates. 

We now study the behavior of these quantities and describe some of their consequences 
in different ranges of QD radii for which aR ^ 1 or aR ^ 1, (corresponding respectively 
to the strong and the weak confinement regimes) and for which we can analytically and 
explicitly compute numerical values. 

III. STRONG CONFINEMENT REGIME 

In this regime, where aR <^ 1, the assumed form of the variational wave function ([3]) 
shall be used and it is appropriate to perform Taylor expansions in Eqs. (jlj) and (jSj) in the 
neighborhood of the dimensionless parameter aR = 0. Because of the analyticity of the 
Taylor expansion of the functions exp(x) and sinh(x), this expansion will remain valid up 
to aR < 1. Thus, we obtain 





where the constants A, B and C are given in closed form by 



2Si 27r - Si 47r 

A = 2 ^ — '- ^ — '- ^ 1.786, 

2tx 




C 



1 



1 



27r2 



0.283. 
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Here, Si(a;) denotes the 



sine integral 
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Consequently, we can deduce an expression of the mean value of the total Hamiltonian H in 
the strong confinement regime as an expansion in powers of a', which is the dimensionless 
variational parameter first introduced in Eq. 

(mid)) .e" _ .„ , E* 



A— - 2{AB - l)E*a' + — a" + ... 



(010) 2/xi?2 kR ^ ' 4 

R 

where the correction terms ". . . " go to zero at least as fast as oc — . The variational 

a* 

parameter a' is now determined by minimizing the expectation value of the energy. Thus, 
we find (Tq = 4:{AB — 1) ^ 0.996 and the corresponding energy value is 

^strong _ ^ ^ _ _ 

2fiR^ kR ^ ' 





-2,R^ 1.T86--0.248E^ (8) 
In the strong confinement regime, this formula has been already analytically obtained |l4l.l33l| 
with trial functions showing the same global form as the previously presented one ([3]) but 
with an interactive part chosen, instead of 0rci(i"eh), equal to 

~ a 
0rcl(rch) = 1 - -'^eh- 

This choice corresponds in fact to the two first terms of the Taylor expansion of Eq. ([3]) in 
the neighborhood of |rch < oR <^ 1. In the meantime, we have succeeded in finding an 
approximate value of its upper bound for R < a*. This bound, as we already mentioned in 
Subsec. Ill CI can be in turn numerically extend to the commonly accepted region of validity 
of QD radii, i.e. i? < 2a*. 



IV. WEAK CONFINEMENT REGIME 



As previously explained, in this regime aR ^ 1, we can retain the global form of the trial 
function used in the strong confinement regime because of the presence of the confinement 
potential. But, as said before, the physical situation requires a modification, before applying 
the variational procedure. 



A. Considerations on the trial function 



Here, the electron-hole pair states consist in binding exciton states, which globally act like 
quasi-particles of total mass M. As matter of fact, in a translation invariant space region, 
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the quasi-particle point of view emphasizes the possibihty for the exciton to show a global 
translational motion in terms of its center-of-mass coordinates. In the weak confinement 
regime, the exciton typical size oc a* is sufficiently smaller than the QD typical radius 
R. So, we can reasonably assume that excitons should exhibit a quasi-particle behavior 
and possess a global translational motion. In order to account for a partial restoration of 
translation invariance, in the weak confinement regime, we have to introduce a center-of- 
mass coordinates plane wave, because the physical behavior of the confined exciton should 
show a reminiscence of the free exciton one, i. e. the wave function of the ground state of the 
exciton shall be considered as the wave function of quasi-particle describing the electron-hole 
pair. Thus, it should have the form 

i){Y^, Th) = ^oio(re)^oio(rh)0rci(rch)0c(rc)0h(rh), (9) 

where 0e,h(i"c,h) = c'^'^°''^ '''=''^. The vectors cr^h are dimensionless quantities which facilitate 
the computation in practice but are not directly physically interpretable. In fact, in order to 
understand properly the meaning of these parameters, we have to define the wave numbers 
vectors respectively in the center-of-mass and in the relative coordinates as 

mlcTc - mlcTh 
'^^^ M • 

Since cTe ■ To + cTh ■ rh = (Tq • fg + (Teh ■ Teh, the functions </'e,h(rc,h) should contribute to 
the exciton total energy by an additional kinetic term, corresponding to the fundamental 
energy of a plane wave in a space region of size R in the center-of-mass coordinates fq, 

1^ |2 \„ |2 \„ 12 |_ 12 _2 

OC 1 = 1 , which IS assumed to be of the form -. This 

ml ml M fi 2MR^ 

constraint is consistent with the conditions 

P=l j 



2_ 1 



(1 + AP 

T72* 

where A = — -. The previous expressions reinforce the cogency of the trial wave function 
m* 

ijj{re, Th). Taking cr^h = does not provide an extra kinetic energy term to the total exciton 
energy in the relative coordinates, except the ones due to the electron and hole confinement. 
But taking IctgP = 1 adds to the electron- hole pair total energy the correct energetic 
contribution in the center-of-mass coordinates, corresponding to the fundamental energy of 
a confined particle of mass M. 
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B. Variational principle 



The trial function ([9]) leaves unchanged the exciton density of probability or the Coulomb 
potential matrix element 

(^AI^A) = (010) and (^/;|l^c(reh)|^) = (0|V^c(reh)|0), 

whereas the free Hamiltonian matrix element gets further contributions, i. e. 

■a 



7r2 E 



2/ii?2 2MR^ 4 

In order to evaluate the mean value of the Coulomb potential in the quantum state 
defined by the wave function tp, we have to compute the following double integral in the 
weak confinement regime for aR > 27r. Since for such QD radii, the convergence of all 
occurring series is insured, we can write 

JJ sin^iTTx) sm'^iny) smh{aRx)e-''^y 

ff dxdy . 2 / TT \ . 2/ ^ \ . , , . _y 

= // sm — — xlsm — —y]smh[X)e ^ 

JJo<x<y<aR ^ v \aR J \aR J 

2am^ ^ k+l[ a^y Jo y ^'"^ r ^ n\ ^ n\ 

27r2 ^ (-47r2)" / An'' \ 



X 

Uo y 

where oF the usual Gauss hypergeometric function. In Eq. 

n>0 ^ 

(ITOl) . we take terms up to the fourth order in the variable in the neighborhood of 0. 

aR 

Therefore, for the terms in Eq. flTUl) whose integrand is not exponentially decreasing for 
y aR ^ 1, we get 



2aRjn y \aRV ^ k + l\ a'^R'J ^ n\ 

-'^ ^ k>0 ^ ^ n=0 
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with 

^' = |{s,(2.)-5iM}«L056, B' = 0, 

C'^-y, D'^^ and E'^\[a-'^^ 

Furthermore, for terms in Eq. flTU]) whose integrand is exponentially decreasing for y — >• 
aR ^ 1, we apply Laplace method (c/. Appendix [A]) to get an answer in the limit of 
aR — > oo. Hence, 



^ • (12) 



In the QD radii range, where the third order actually contributes to Eq. f|TT]) . i.e. for radii 
near the lower bound of the possible radii in the weak confinement regime aR ^ 27r — 
which, as we will see later, corresponds to a size R ~ vra* — , the third and fourth order 
terms contribute to the exciton energy with terms of the same order of magnitude. But for 

n 

radii larger than 4a* p^, the third and also the fourth order terms are in fact irrelevant. 
So, we can stop the expansion in Eq. f|TT]) at the second order term. Thus, the Coulomb 
interaction diagonal matrix element and the square of the norm become 



^ 8A' ( 4 7r2 D" tt^ 



kR a2i?2 3 a^R2 A' a^R^ \ a^R"^ 



a^R^ 3 (T2i?2 2 A' a^R^ \ a^R'^ 



(13) 



E' 

where by definition D" = ~ 0.372. Finally, we obtain in this regime the following 
cr'-expansion of the variational energy 

mm 7r2 E* ,2 , 2 1 



+ —a'^ - E*a 



{^ip) 2/ii?2 2Mi?2 4 3/ii?2a' 
3/}" 7r2 1 a* 

H TT h . . . 

4 A' i^R^a'^R 

Applying now the variational principle, we get the value of the dimensionless variational 
parameter a' 

2 ofa*\^ 3D" ,/a*V 
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Because the second and the third terms in the previous expression do not contribute to the 
total electron-hole pair energy, we conclude that ctq ~ 2 and 



= -E*-\ \ (lA) 

6fiR^ 2M{R - r/(A)a*)2 ' ^ ^ 

from which we are able to extract an explicit expression for the Kayanuma-function as 

„(A) = 6 ^ ' with 5 = — — TT = — 1 ]^ 0.208 and A = As we have 

^ A 16 A' 40 V 8A' ) ml 

already specified, let us observe that ctq ~ 2 implies that R > 7ra*, which is surprisingly 
consistent with the usually accepted numerical region of validity. Moreover, we see that the 
interaction part f/'rcilrch) of the wave function coincides with the ground state wave function 
of a reduced mass yU hydrogen-like atom up to a normalization factor and contributes to 
the excitonic total energy with the leading binding fundamental energetic term —E*^ as we 
expect. 



The ?7-function was first introduced phenomenologically by Eq. (28) in jlj] to get a better 
fit between numerical and empirical results; but no analytic derivation of it is available so far. 
First, here, we manage to give an expression for this function, which satisfies the electron- 
hole exchange symmetry: A — > A~^. Second, as was already depicted by Kayanuma, the 



QD size renormalization term ri{X)a* is interpretable as the so-called dead layer 3J]: this 
is the physical reminiscence of the fact that, although it could be successfully described as 
a quasi-particle, the exciton is not actually itself an indivisible particle. In fact, its center- 
of-mass motion, the one on which the quantization is really performed, could not reach 
the infinite potential well surface unless the exciton undergoes a strong deformation in the 
relative motion of the electron and the hole. This implies that the picture of a point-like 
exciton is no longer appropriate in this region of space. The most convenient way to study 
this exciton consists in thinking of it as a rigid sphere of radius oc a*, where the proportional 
factor must not exceed | too much, because /* = |a* is the mean value of the relative 
distance between the electron and the hole in the non-confined exciton ground state. We 
have to get the largest possible radius I* for the sphere picture, when A oo, i.e. in 
the limit of infinite hole mass because the hole stays motionless at the center-of-mass of 
the electron-hole system. Furthermore, the smallest possible radius must be obtained in 
the symmetrical case A = 1 and must be held at |/* = |a*, which imposes 77(1) ~ 0.75. 
This matches quite well with both experimental and theoretical results as shown in Tab. [B 
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Finally, regarding the different values of 77(A) around | at A = 5, we can safely conclude 
that the picture of infinite hole mass becomes valid as soon as A ~ 5. 



Table I: Comparison between numerical and analytical values for the r^-function. 



A 


1 


3 


5 


f?num(A) 


0.73 


1.1 


1.4 


??theo(A) 


0.83 


1.1 


1.5 


relative error 


14% 


<1% 


7% 



In Tab. [H we compare numerical values computed from numerical simulations taken by 
this function 77(A) for A = 1,3,5 to those theoretically predicted by the previous relation 
and note that there exists a reasonable agreement between both results. We even succeed 
to enlarge the region, in which the weak confinement regime is satisfactory. 

V. A PSEUDO-POTENTIAL-LIKE METHOD 

The problem, that we still have to deal with, consists in finding a physical way to substract 
off the term reminiscent of the kinetic energy term in Eq. f|T^ . As matter of fact, because of 
the presence of the reduced mass, it is interpretable as a kinetic energy term in the relative 
coordinates. However, this type of exciton kinetic energy has to be already contained in the 
Rydberg energy term, because of the validity of the Virial theorem in the relative coordinates. 
Besides, in this assumption, the higher order contributions to the electron-hole energy must 
be interpreted as higher order contributions to the kinetic energy of the exciton, viewed as 
a quasi-particle of mass M, which physically justifies the intuitive idea that for very large 
radius R only the quasi-particle point of view should be responsible for the exciton kinetic 
energy. 

A. Weak confinement regime 

In analogy with the pseudo-potential method for metals, we propose to introduce an 
additional potential term W in the Hamiltonian H, which shall depend only on the electron- 



15 



hole relative distance r^h, so that it contributes to the second order of the exciton total energy 
but not to the third order — i. e. the one which allows us to determine the expectation value 
of the ?7-function. Furthermore, we assume that W must 

• be attractive at distances in the range of a* in order to promote excitonic state with 
typical size around his Bohr radius; 

• be repulsive at short distances in order to penalize excitonic state with small size and 
to remind that the exciton typical size should not be smaller that the typical size of 
the surrounding lattice spacing; 

• be exponentially small for large distances, in order to not perturb the long range 
Coulomb interaction. 

Consequently a natural choice of such potential should be of the form 



Inspection shows that it satisfies all of the previous constraints; and, in the weak confinement 
regime, it exactly provides us what we ask for 



Remark. The potential amplitude in Eq. f|T5|) is fixed ad hoc in order to get the previ- 
ous correct kinetic energy contribution. Consequently, the pseudo-potential form we choose 
seems to be arbitrary. However, the form of W{rch) oc ?"eh^~^"^ only one which 

straightforwardly first contributes to the second order term in both weak and strong con- 
finement regimes and which does not change the exciton energy behavior to zeroth and first 
order terms. 

B. Strong confinement regime 

As we have introduced a new term in the exciton total Hamiltonian H, we have to check 
its consequences in the strong confinement regime. The most important one is the significant 
decrease of the expected value of the exciton energy because of the exciton total Hamiltonian 



327r2 
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redefinition: H' 

as 



= H + W. In fact, the pseudo-potential W contributes to the second order 



(0|iy(reh)|0) 
(010) 



CEU 1 



'~9~ 



-19.9E* + ... 



f(4 



This expression is deduced by a reasoning similar to that of Sec. IIIII Therefore, we must 
keep in mind that it is only valid when 2R < a* because of the pseudo-potential exponential 
dependence, which changes the validity conditions of the Taylor expansion we made. 

Figure 1: Behavior of the excitonic ground state energy as a function of the QD radius computed for 
a confining infinite potential well with ( — ) or without ( — ) the presence of the pseudo-potential W 



and for a confining finite potential step of height Vq ~ leV ( ) 



results for CdS 



21[ and compared to experimental 



35[ microcrystallites with material parameters: e = 11.6, m* = 0.235mc, 



1.35me, Eg = 2.583eV, E* = 27meV and a* = SO.lA, where rric is the electron bare mass. 
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As we can see from Fig. [T], the excitonic energy computed in the presence of the pseudo- 
potential gets a significant hetteT fit to experimental results, in the validity domain 2R < a*, 
than the one calculated without this tool. Nevertheless, the divergence for very small QD 
size still holds as a consequence of the infinite potential well assumption. In order to extend 
the validity domain upper bound, it would be efficient to carry on the energy expansion 
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to one or maybe to two further orders. But, at this point, the computation become quite 
involved and we cannot insure the success or even the relevance of this approach. 



VI. CONCLUSION 

By using an improved EMA model, we are able to obtain well known results in the 
strong confinement regime, to correctly apprehend the weak confinement regime by adding 
a pseudo-potential term to the total Hamiltonian describing the electron-hole dynamics and 
to find an nice analytical expression for the phenomenological function t] of Kayanuma. 

As we have shown, the introduction of the pseudo-potential W allows to reduce in part 



the overestimation due to the infinite potential well observed in [2l|, particularly near the 
upper boundary 2R ^ a*, even if it still does not accurately approximate the total energy 
for much smaller QDs. For such QD sizes, the best approach would probably require a finite 
confinement potential. Thus, a future research work might focus on applying this pseudo- 
potential method, in the presence of a confining finite potential step, in order to obtain a 
hopefully better behavior for excitons in very small QDs. 

Appendix A: THE LAPLACE METHOD 

The purpose of this method is to study asymptotic behavior of integrals like I{t) = 

/ dxg{x)e*^^^^ and rigorously determine a mathematical equivalent for such quantities when 
Jo 

t ^ oo. 



One can prove the following theorem 36|: Let a > and g, h : [0, a] — > M two continuous 
mappings such that 

i. / dx \g{x)\e'^^''^ < oo; 
Jo 

ii. 3(5o > / V(5 G [0, 6o] and Va; G [6, a[, h{x) < h{6y, 
iii. g{x)^ Ax°' and h{x) = b — cx^ + o(x^), where a > — 1 and c, /3 > 0. 

a:0+ 
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Then, I{t)^-T 



too 



( 



a + 1 
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